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PREFACE 



\ 

; 



This study was initiated as Subtask 1, TNT Equivalency Study to NASA Study 
C-II, Advanced Missions Safety Studies. Other studies in this series are 
Subtask 2, Safety Analysis of Parallel versus Series Propellant Loading of the 
Space Shuttle, Aerospace Report No. ATR-7 1(7233)- 1 and Subtask 3, Orbiting 
Propellant Depot Safety Study, Aerospace Report No. ATR-71(7233)-3. 

This study was supported by NASA Headquarters and managed by the Advanced 
Missions Office of the Office of Manned Space Flight. Mr. Herbert Schaefer, 
the Study Monitor, supported by Mr. Charles W. Childs of the NASA Safety 
Office, provided guidance and counsel that significantly aided this effort. 

Study results are presented in three volumes; these volumes are summarized 
as follows: 

Volume I: Management Summary Report presents a brief, concise 
review of the study content and summarizes the principal conclusions 
and recommendations. 

Volume II: Technical Discussion provides a discussion of the 

available test dat? and the data analysis. Details of an analysis of 
possible vehicle static failure modes and an assessment of their 
explosive potentials are included. Design and procedural criteria 
are suggested to minimize the occurrence of an explosive failure. 

Volume III: Appendices contains supporting analyses and backup 
material. 
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APPENDIX A 


Tabulation of Pyro Test Data with 
Yield in Percent TNT as Determined 
by URS Systems Corp. , Bellcomm Inc . , and The Aerospace Corp. 



In this appendix, the following table headings are defined as: 

L/D t - Tank Length /Diameter 

D /D - Rupture Diaphragm Diameter /Tank Diameter 
o t 

D/W*^ - Distance from Event/(Total Wt of Propellants)*^ 

The legend below is a key for the Ignition Source column. 

A Bottom tank did not release fluid until after impact of top tank 
B Tank rupture ignition 
C Cap ignition at tank rupture 
D Cap ignition 
E Diaphragm rupture ignition 
F Cap ignition at 100 psi internal pressure 
G Ignition caused by poor venting 
H Fire on tank before ignition 
I Squib ignition 
J Probable impact ignition 

K Fire visible at top of tank 10-15 ms before ignition 
L. Ignition at impact 
M Fire at top of tank 

N Self-ignition - suspect diaphragm break 
P Fire at top of tank after rupture, ignition at bottom of tank 
Q Tank fell through stopper 
R Fire on pad before ignition 
S Tank broke apart at impact 
T Self- ignition - tank leaked 

U Cap ignited after propellant on pad. Fire went out; detonation at 
time given. 
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Summary of Test Methods Employed by Investigators 
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GENERAL 


B. 1 

The data produced by these tests are summarized and analyzed in Volume II. 

The purpose of this appendix is to present the configurations and pertinent 
parameters used in the test programs. 

B. 2 A. D. LITTLE, INC., SPILL TESTS (Ref. B-l) 

A. D. Little, Inc. conducted a series of LO^/LH^ spill tests in 1962. All 
tests were conducted with an oxidizer-to-fuel ratio of 5 to 1. Nine tests in- 
volved a total propellant weight of 45 lb each, and three of 225 lb each. 

As shown in Figure B-l, tanks containing the LO^ and LH^ were tipped toward 
each other, spilling the propellants in partially impinging streams into a 
shallow pit to achieve active mixing. Tank release was sequenced to ensure 
that both tanks were exhausted simultaneously. An explosive charge in the 
spill pit was used to initiate the explosion in all tests. 

Overpressure was sensed by piezoelectric pressure transducers, and was 
recorded by photographing an oscillograph trace. Five gauges were positioned 
along a single radial line from the explosion. Distances varied from 30 to 
70 ft in the 45-lb tests, and from 60 to 120 ft in the 225-lb tests. 

Test results were reported in terms of pressure yield, calculated by averaging 
the pressure yields determined individually for each data point. The TNT 
reference curve used in determining yield was developed experimentally by 

A. D. Little, since their test program preceded the publication of the currently 
accepted TNT curves prepared by Ballistic Research Labs. (BRL) (Ref. B-2). 

B. 3 NASA SPILL AND TANK TESTS (Ref. B-3) 

In 1964, NASA-MFSC conducted a series of LO^/LHj, Each of these 

tests involved 200 lb of propellants, and an oxidizer-to-fuel ratio of 5 to 1. 

Six spill tests and seven tank tests were run. 
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The description of the test method employed in the A. D. Little spill test is 
applicable also to the NASA spill test method (Figure B-l). In the tank tests, 
the propellants were contained in tandem tanks, separated by an intertank 
bulkhead (Figure B-2). The simulated failure modes included intertank bulk- 
head rupture and tank wall rupture. Bulkhead rupture was produced either 
by overpressurization or by detonating a linear- shaped charge attached to the 
structure. Tank wall rupture was accomplished by a linear- shaped charge. 

It was intended to control ignition delay by utilizing either a squib (for flame 
ignition) or a blasting cap (for shock ignition). However, in all tests ignition 
was spontaneous at tank failure and was probably caused by the shock accom- 
panying the tank failure. 

Overpressure was sensed by piezoelectric pressure transducers and recorded 
on magnetic tape. Six gauges were positioned along a radral line at distances 
ranging from 25 to 250 ft. Three additional gauges were located in a radial 
line 45 deg away from the first line, and at distances of 25 to 65 ft from the 
test article. Four gauges were spaced along a third radial line, 80 deg away 
from the first line, at distances of 25 to 100 ft. 

For each series of tests a range of pressure and impulse yields was reported. 
The ranges were determined by plotting the high and low TNT yield percentage 
curves which envelop all of the data points. The BRL-TNT curves were used 
as the reference basis. No tabulated test data were reported. 

B. 4 AEROJET-GENERAL CORP. CONTROLLED INTERFACE 

TESTS (Ref. B-4) 

In 1965, AGC conducted a series of LO^/LH^ tests to evaluate the explosive 
characteristics of these propellants. The propellants were combined using a 
mixing technique that permitted controlled variation of the oxidizer/fuel 
contact area and the weight of propellant. A fixed 5 to 1 oxidizer-to-fuel 
ratio was used for all of these tests. The propellants were mixed by placing 
oxidizer -filled glass dewars in a fuel-filled pan and shattering the dewars 
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with an explosively generated shockwave (Figure B-3). Two tests were run 

at each of three weights: 100, 150, and 2Z5 lb. A constant 4 to 1 ratio of 

? 

contact area between ‘ he oxidizer and fi el (in ft ) to total propellant weight 
(in lb) was maintained for all tests. 

Overpressure was sensed by piezoelectric transducers and recorded on 
magnetic tape. Four gauges were positioned along each of three radial lines 
spaced 120 deg apart. Gauge distances from the test article varied from 25 
to 98 ft. 

Test result* were reported in terms of pressure yield. Yield was determined 
by superimposing the test data points on plot;, of TNT percentage yields. T>e 
TNT reference curves were determined experimentally and were reported as 
being in good agreement with the BRL curves. Tabulated test data were 
reported. 

B. 5 URS SYSTEMS CORP. PROJECT PYRO TESTS (Ref. B-5) 

) 

Project Pyro, conducted in the 1965-1967 time period, had as its objective 
the development of a reliable philosophy for predicting the credible damage 
potential which may be experienced from the accidental explosion of liquid 
propellants. Two general failure modes were simulated. In the Confinement 
by Missile mode (CBM) failure occurs in the intertank bulkhead, and all 
propellant mixing occurs within the tanks. The Confinement by Ground Surface 
(CBGS) condition simulated the case where the propellants spill from the tanks 
and mix on the ground surface. Both vertical (CBGS-V) and horizontal 
(CBGS-K) spill cases were simulated. A total of 78 tests were conducted and 
raw data were reported; 28 were in the CBM mode, 39 were CBGS-V and 1 1 
were CBGS-H. The variations in test article parameters and configurations 
for the three modes, as well as the general arrangements of the test setups, 
are depicted in Figures B-4, B-5, and B-6. 

Overpressure sensing was accomplished with piezoelectric transducers, and 
data were recorded on magnetic tape. Three radial lines of gauges were used, 

1 
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• LI^ in Aluminum Container 

• LO£ in 1 Liter Glass Dewars 

• Oxidizer/ Fuel Weight Ratio 
5: 1 (Constant) 

2 

• 0. 25 ft Contact Area per 
Pound of Propellant (Constant) 

• Dewars Shattered by Mild 
Detonating Fuse Under the 
Container Bottom 



No. of 
Tests 

Total 

Propellant 
Wt. (lb) 

Container 
Size (in. ) 

No. of 
Dewars 

Contact 
Area (ft^) 

2 

100 

12 X 24 X 24 

36 

25. 1 

2 

150 

12 X V- Jb 

52 

36. 8 

2 

225 

12 X 44 X 44 

81 

56. 2 


Figure B-3. Aerojet-General Corporation Test Program 
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^ ( 

; spaced 120 deg apart. Five transducers were positioned along each line. 

Distances varied from 23 to 67 ft for the 200 and 1, 000-lb tests, and from 
38 to 600 ft for the 25, 000 and 91 , 000-lb tests. 


> Test results were reported as terminal yield, which was apparently an 

? average of the pressure and positive phase impulse yields. The BRL-TNT 

\ pressure and impulse curves were used as a reference. Tabulated data for 

I each test were reported. 
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Figure B-4. Project Pyro Test Program, Confined by Missile 




Tandem tanks with 5 mil aluminum foil 
bottoms 

Tanks dropped to within 1-1/2 tank 
diameters of pad 

Tank bottoms ruptured simultaneously 
by star cutters 

Propellants spill onto pad, forming 
overlapping pools 

Planned ignition by linear charge 

LO^ / L.H2 weight ratio 5:1 

Parameters varied: 

Total propellant wt. (200, 1000, and 
25, 000 lb) 

L/D ratio (1.8 and 5.0: 1) 

Impact velocity (23, 44, and 78 fps) 


Figure B-5. Project Pyro Test Program, Confined by Ground 
Surface - V ertical 


B-* 



• Separate tanks with aluminum 
foil bottoms 

• Tanks dropped separately - bottom 
tank stopped within 1-1/2 diameters 
of pad 

• Tank bottoms ruptured by star 
cutters 

® Propellants spill onto pad, forming 
overlapping pools 

• Planned ignition by linear charge 

• LO z / LH 2 weight ratio - 5:1 

• All tests 200 lb total propellant 
weight 

• LO 2 tank impact velocity 12 fps 
(all tests) 

• LH 2 tank impact velocities 23 fps 
and 78 fps 

• L/D ratio 1.8:1 (all tests) 


Figure B-6. Project Pyro Test Program, Confined by 
Ground Surface -Horizontal 
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GENERAL 


C. 1 

This appendix contains a brief synopsis of studies conducted at the University 
of Florida under the direction of Dr. E. A. Farber. 

The studies to be discussed are: the construction of a mathematical model 
for predicting explosive yields of liquid propellants, a seven chart approach 
for systematically applying the model, and the critical mass phenomena. 

These studies are thoroughly documented in Refs. C-l through C-5. 

C. 2 MATHEMATICAL MODEL 

The mathematical model begins with the assumption that the explosive yield 
function is related to a mixing function where the yield function, y, is defined 
as a fraction of maximum theoretical yield, and the mixing function, x, is 
defined as the fraction of propellants actually mixed. By dealing with the 
yield function and the mixing function in this way, a generalized expression can 
be developed which is assumed to be applicable to all liquid propellants. 

Using test data produced by A. D. Little (Ref. C-6), Figure C-l was constructed 
and a least squares fit of the data points shown. Data points Jl, J2, J3, are 
from tests made with a l/25th scale Saturn V vehicle, while data points D and 
H were maximum values observed in smaller scale testing of a Saturn V 
configuration. In each case, the quantity of propellants mixed at the time of 
ignition was known so that the fraction of propcTlants mixed could be plotted 
against the fraction of the theoretical yield that was observed. The plot of 
the data suggests a relationship between the fraction of theoretical yield 
and the fraction of propellants mixed of the form y = C(x)^ where 0sx5 1, 

0 < y < k, and k < 1 . The yield function and the mixing function each are 
presumed to be a random variable. 
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Moving forward with the assumption that the mathematical model does 
represent a generalized statement of explosive yield, a mathematical problem 
was addressed to find a bi-variant probability density function where the 
expected mean value of the yield for any given mixing function value is: 

E(x, y) = C* 


o 


s 

i o 


■ 

•v 


0 f(x, y)dy 


The following bi -variant function was selected for x and y. 


.. . dr(a + b + c) d-1,. d.a-1 b-1. d »c-l 

f(x ’ y> ° x (1_x) y (x - y) 


The only restrictions on this equation are given as y > 0, x > 0, y 5 x^, and 
d * 0. 

Mathematical relations were defined which permitted estimating the values of 
parameters a, b, c and d using experimental explosion data. It was determined 
that three of these four controlling parameters are constants: b = 4. 0, c = 1.1, 
d = 1.5. Parameter a is variable and it was demonstrated that it varies with 
the propellant weight and is thus a scaling function applicable to the model. 

When the controlling parameters of the probability distribution have been 
estimated, it is possible to determine the probability distribution for the yield 
function, the probability distribution for the mixing function, the combined 
yield-mixing probability distribution, and confidence limits for the yield 
function, and the mixing function can be assigned. The typical shape of the 
mathematical model represented by a statistical surface is illustrated in 
Figure C-2 and is likened to a shark fin. 
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Figure C-2. Statistical Surface of Explosive Yield Model 
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The relationships of the scaling factor, parameter a, to propellant weight 
and to explosive yield are shown in Figure C-3. The smaller the value of 
parameter a, the higher the yield, and, based on existing data points, higher 
propellant weights are associated with smaller values of parameter a. It 
must be emphasized that these two relationships have been developed using 
experimental explosion and accidental missile failure data in order to estimate 
the value of parameter a. The model is so completely flexible that oti er values 
of parameter a could derive from different experimental or missile fail ire data. 
Different values of a could substantially change predicted yields. The greater 
the size of the data bank, the greater the confidence that may be associated 
with yield predictions from this mathematical model. 


A plot of explosive yield as a fraction of the potential yield versus propellant 
weight is shown in Figure C-4; also the 95% probability limit is indicated. 
Superimposed are data plotted from Project Pyro and from a number >f actual 
missile failures shown in Table C-l, where the explosive yield has usually been 
estimated on the basis of analysis of resulting damage. It can be said that 
the available data fit we)' within the prediction limits presented. It is also 
evident that the fit is obtained because these data were used to determine the 
controlling parameters of the prediction equation. Again, the lack of data 
points above approximately 250, 000 lb of propellants gives rise to some doubt 
as to the propriety of extrapolating to propellant weights in the millions of 
pounds range. 

C. 3 SfcVSN CHART APPROACH 

A systematic procedure for the analytical prediction of the yield from liquid 
propellants was also developed. The procedure is summarized in seven charts 
(Refs. C-l, C-4, and C-5). This method has the potential to produce more 
information concerning explosions and their yield, but it also requires 
considerably more input information and knowledge about the liquid propellants 
and the mechanisms contributing to the explosion. The input information is 
not presently available in many cases and therefore was assumed in developing 
the theory. 
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Table C-l. Yield Estimates and Data of Missile Failures 
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In the Seven Chart Approach, the complex question of explosion yield 
prediction was divided into three parts which were examined separately and 
independently and then combined to give the desired results. These parts are 
discussed in the following paragraphs. 

C. 3. 1 The Yield Potential 

The yield potential is defined as the maximum theoretical explosive yield which 
can be produced at any time t if the propellants present are mixed in the most 
favorable manner. This function can be calculated for any given propellants 
and selected oxidizer -to -fuel ratio as a function of time using chemical 
kinetics, heat transfer theory, knowledge of the failure mode and the original 
configuration of the propellant containers. 

C. 3. 2 Mixing Function 

The mixing function is defined as the fraction of the propellants which is actually 
mixed at any time t. It is only this mixed portion of the propellants which 
takes part in an explosive yield. The fraction can be defined in terms of 
volume modified by factors which define the degree of mixing which has taken 
place, or in terms of contact area between the propellants and/or other like 
considerations. The mixing function is controlled by the type of propellants, 
the container size and configuration, and the failure mode. The mixing function 
always increases from zero at the time of failure to a maximum and then 
decreases again if ignition does not occur. The mixing function multiplied by 
the yield function gives the expected yield as a function of time for the 
propellants, configuration, and failure mode under consideration. 

C. 3. 3 Ignition Time 

The ignition time determines what the predicted explosive yield will be since 
it defines the specific point on the expected yield versus time curves generated. 
Early ignition results in low yields as does very late ignition. Somewhere 
between, optimum conditions occur which will produce the maximum explosive 
yield prediction for the conditions being considered. In the case of the actual 
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vehicle failures which involve large quantities of propellants, there is the 
possibility of having early ignition, since many ignition sources would exist, 
e. g. , electrostatic charges, failures of electrical systems, failure of 
structural members, hot engine parts, etc. , which would result in low 
explosive yields. 

C. 3. 4 Conclusions 

It is obvious that the calculation of the yield potential is a complex exercise. 
Similarly, arriving at the mixing function is equally as difficult. Laboratory 
experiments have been conducted which were aimed at verifying necessary 
assumptions. Pellets have been vibration -mixed, wax casts have been made 
of the mixing process, high speed photographic records have been made of the 
mixing process, and high sensitivity thermocouple grid records during the 
mixing period (up to de struct by the explosion) have been analyzed. It is 
suggested by Dr. Farber that when experimental explosion results are not 
available or feasible to obtain, the Seven Chart Approach will provide a 
prediction technique which takes into account the kinds of propellants, vehicle 
configuration, and failure mode. While the experiments tend to support the 
theory, not enough data are available yet from use of the actual propellants in 
simulated scale configuration or in sufficiently large quantities to be conclusive. 

C. 4 CRITICAL MASS PHENOMENA 

The Seven Chart Approach indicates that the explosive yield is a function of 
the ignition time, which is that period from the start of mixing until ignition 
occurs. Dr. Farber and his associates were impressed that early ignition 
occurred frequently with large liquid propellant quantities. Although many 
sources of ignition are available during a vehicle failure, they considered 
the possibility that the mixing processes themselves produced ignition as the 
result of electrostatic charges generated and discharged across vapor bubbles. 
Laboratory tests have been conducted confirming in fact that an electrical field 
is generated during the mixing process. Calculations have indicated that a 
critical mass of about 2300 lb of LO^/LI^ would produce a sufficiently strong 
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charge capable of sparking, thereby causing ignition to occur. These 
calculations presumed that the mixing was due primarily to boiling of the 
propellants. If the propellants were brought together more violently, more 
propellants could be mixed before the voltage build-up occurred. Similarly, 
if the propellants were mixed more gently, greater quantities could be mixed. 
It is expected that additional work on the theory of critical mass will be 
forthcoming in the near future. 
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GENERAL 


D. 1 

The current configuration of the space shuttle in fully loaded condition on the 
pad contains approximately 4 million lb of LH 2 and LC> 2 propellants. Of these, 
2. 9 million lb of L0 2 and 0. 5 million lb of LH 2 are stored in two tanks aboard 
the booster. The orbiter has two L0 2 tanks with a capacity of 0. 25 million lb 
each and a single LH 2 tank with 0. 1 million lb capacity. 

The current criteria specify that the TNT equivalence for L0 2 >LH 2 propellant 
combinations is 60%. Based on this equivalency, the facilities siting for the 
space shuttle would have to be based on a total yield equivalent to 2. 4 million 
lb of TNT. 

Studies for other vehicles have indicated that the probability of a failure 
occurring, such as tank rupture, which could result in the explosion of 
propellants such as LH 2 and L0 2> is very low. For instance, data from an 
analysis for the Titan HIM vehicle indicated that the probability of a tank 
rupture during the time period from T-30 to T-0 minutes was approximately 
0. 4 X 10' 6 . A specific failure mode and effects analysis has not been made 
for the space shuttle vehicle; however, it is generally considered that the 
probability of a tank rupture while on the pad would also be very low. 

In spite of the anticipated low probability of a situation occurring which could 
result in an explosion on the pad, there is the question of the magnitude of the 
explosion if such an event occurs. This appendix presents a statistical 
development of the yield of exploding propellants for a multi -tank vehicle 
configuration such as is used on the space shuttle. 

In that analysis, a probability density function for the yield of an explosion was 
first developed. The properties of that probability density function were then 
used to establish a probability density function considering the fact that various 
combinations of tanks of the space shuttle (hence, different quantities of 
propellants) may be involved in the explosion. The composite probability 
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density function then was used as the basis for probability statements 
relative to the /ield of the space shuttle propellants. 


3 


D. 2 APPROACH 

In the event of an explosion, the yield, y, may vary from zero to some 
maximum value, y^. Many factors obviously affect the yield of the explosion, 
including the mixing of the propellants. However, adequate information was 
not available in this study from which to base the development of a probability 
density function for yield. Therefore, several assumed probability density 
functions were evaluated in this analysis to show the sensitivity of the results 
to this input. The various density functions used are discussed in subsequent 
sections of this appendix. 

A basic requirement for each of the probability density functions evaluated was 
that the following mathematical condition be satisfied: 

1 

F(z)dz = 1 

where f(z) is the probability density function whose value must be equal to or 
greater than zero for all admissible values of z and z = 7/y^. 

Using the assumed distributions on yield, a composite probability density 
function was then developed as discussed below. 

The ST 3 vehicle contains five tanks for which there are 31 possible combinations 
which may be involved in an explosion. To establish the probability density 
function tor the complete system of tanks, the following ground rules were used: 

a. When the ruptured tank(s) is an tank, it was assumed that 
an explosion will x?.ot occur but that a fire may result. 
Subsequent explosions of the other propellants aboard th* 
vehicle as a result of the fire will be low order and will not 
be c nsidered. 
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b. In order to have an explosion involving all the propellants 
of two or more tanks (one or more of which must contain 
LiH 2)> the ruptures of the individual tanks have to occur 
essentially simultaneously to result in a high order explosion 
rather than a fire. 

c. A rupture of a tank induced by impact of debris resulting from 
the rupture or explosion of another tank would result in fire or 
subsequent lower order explosions than those associated with 
the initial rupture and therefore need not be considered 

The above assumptions result in a widely varying range of quantities of propel- 
lant which may be involved in the explosion. The various combinations of tanks 
are shown in Table D-l together with the amount of propellant in each tank and 
the total propellants that would be involved in each case. 

Since there are five tanks in the booster and orbiter, the probability that each 
of the possible combinations of tanks shown in Table D-l will be involved in the 
event must be considered in establishing a probability density function for the 
complete system of five tanks. Using the above ground rules and the assump- 
tion that the tanks are equally likely to fail, the conditional probability that a 
specific tank will be involved at the time of the rupture is 0. 2. The probability 
that any two specific tanks would be involved was assumed to be (0. 2)(0. 2) - 
0. 04, etc. Probability values so obtained were then used with the data in 
Table D-l and the previously discussed distribution functions to establish the 
composite probability density function for the complete system. 

The equation for the composite probability density function is 


f(7) = 


dp„ F n (y)] 


where 

p = the probability that tank combinatior u (as defined in Table D-l) 
n is involved in the initial explosion. 

F n (y) = the density function for tank combination n. 
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Table D-l. Propellant Tank Combin-tions and Capacities (lb X 10^) 


Tank 

Combination, 

n 

Booster 

Orbite r 

Total 

Total 

Propellants 
Contributing 
to Explosion, 


o 

H 

o 

O 

H 

O/H 


i 

2. 9 





2. 9 

0. 0 

2 


0. 5 




0. 5 

0. 5 

3 



0. 25 



0. 25 

0.0 

4 




0. 25 


0. 25 

0. 0 

5 





0. 1 

0. 1 

0. 1 

6 

2. 9 

0. 5 




3. 4 

3. 4 

7 

2. 9 


0. 25 



3. 15 

0. 0 

8 

2.9 



0. 25 


3. 15 

0. 0 

Q 

/ 

2. 9 




0. 1 

3. 0 

3. 0 

10 


0. 5 

0. 25 



0. 75 

0. 75 

li 


0. 5 


0. 25 


0. 75 

0. 75 

12 


0. 5 



0. 1 

0. 60 

0. 60 

13 



0. 25 

0. 25 


0. 50 

0. 0 

14 



0. 25 


0. 1 

0. 35 

0. 35 

15 




0. 25 

0. 1 

0. 35 

0. 35 

16 

2. 9 

0. 5 

0. 25 



3. 65 

3. 65 

17 

2. 9 

0. 5 


0. 25 


3. 65 

3. 65 

18 

2. 9 

0. 5 



0. 1 

3. 00 

3. 50 

19 

2. 9 


0. 25 

0. 25 


3. 40 

0. 0 

20 

2. 9 


0. 25 


0. 1 

3. 25 

3.25 

21 

2. 9 



0.25 

0. 1 

3. 25 

3. 25 

22 


0. 5 

0. 25 

0. 25 


1. 0 

1. 0 

23 


0. 5 

0. 25 


0. 1 

0. 85 

0. 85 

24 


0. 5 


0. 25 

0. 1 

0. 85 

0. 85 

25 



0. 25 

0. 25 

0. 1 

0. 60 

0. 60 

26 

2. 9 

0. 5 

0. 25 

0. 25 


3. 90 

3.90 

27 

2. 9 

0. 5 

0. 25 



3. 75 

3. 75 

28 

2.9 

0. 5 


0. 25 

0. 1 

3. 75 

3.75 

29 

2. 9 


0. 25 

0.25 

0. 1 

3. 50 

3. 50 

30 


0. 5 

0. 25 

0.25 

0. i 

1. 10 

1. 10 

31 

2. 9 

0. 5 

0. 25 

0. 25 

0. 1 

4. 0 

4. 0 
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D. 3 


RESULTS 


) 


D. 3. i Model 1 

As indicated previously, several probability density function models were 
evaluated in this study. For Model 1, the following criteria were used as the 
basis for selection of an appropriate density function: 

a. The value of 60% TNT equivalency was assumed to be the 
mean value, y, of the probability density function. It was 
also assumed that all yields in the neighborhood of 60% 
are highly likely on a relative basis. 

b. All values of F(z) in the vicinity of zero yield and the 
maximum value are relatively low in comparison with the 
value of F(z) In the vicinity of the mode. 

c. It was assumed that the maximum value, y n , would be based 
on a TNT equivalency of 100%. That is, the explosive effects 
would never exceed that for a similar amount of TNT. 

Using the above conditions as criteria, a Beta function was selected which 
has the general equation: 


F(z) = 


(o- + p + 1)! 
al p ! 


z a (l - z) 


P 


where 

z = normalized form of the random variable (yield) 

= y'Vn 

The Beta function describes a family of curves which is dependent upon a 
choice of a and |3. For this analysis, the desired characteristics of the 
distribution, for the case where the maximum yield is 100% eauivalency, were 
obtained using a value of a ~ 5. Using this value and the relationship for the 
mean 

— a + 1 
z = 




* 


See Ref. 


B-l for a detailed discussion of the Beta function. 
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the value of (3 = 3 was determined. The equation for the distribution function 
is then 


F(z) = 504 z 5 (l - z) 3 


A plot of this equation is shown in Figure D-l. This curve indicates that the 
most likely equivalency (modal value) is 0. 625, whereas the mean value was 
0. 60. It also shows that there is relatively low probability of equivalencies 
less ti — n 0. 2. This density function in terms of yield is given by the equation 


F(y) - 



3 


The results of the analysis using the density function shown in Figure D-l 
are presented in Figure D-2. The unusual shape of this composite density 
function at low values of yield results from the basic shape of the probability 
density function for an individual explosion which was assumed and the low 
values of yield which result from many combinations of the tanks. This figure 
shows that the most likely value of yield is less than 0. 1 X 10^ lb while the 
expected yield is 0.56X 10^ lb. It is also interesting to note that the 95% 
cumulative probability occurs at 2. 3 X 10^ lb. In Figure D-3, the probability 
of exceeding a given specific yield is plotted. 

D. 3.2 Model 2 

In Figure D-4, the second model considered is shown. In this case the mean 
value of the distribution of 60% was retained but the maximum value was 
increased to 150% TNT equivalency to acknowledge the possibility of higher 
yields than were possible with Model 1. 

The composite probability density function and the cumulative probability are 
plotted in Figures D-5 and D-6. They show that the expected value remains 
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Figure D-2. Composite Probability Density function - Model 1 





Figure D-3. Probability of Exceeding a Given Yield - Model 1 












approximately the same as for Model 1 while the 95% cumulative value has 
increased slightly due to the fact that higher yields are possible using this 
model. 

D. 3. 3 Model 3 

In Figure D-7, the probability density function is shewn for Model 3. In this 
model, the mean value of the distribution was assumed to be 20% compared 
to 60% in the previous analyses. The other general characteristics of the 
distribution are similar to those in the preceding case. 

In Figures D-8 and D-9, the results of the analysis for this model are 
presented. They show that both the expected value and the 95% cumulative 
value are reduced substantially from the previous models considered. 

D. 3. 4 Model 4 

Available experimental data were also used to develop a density function. In 
Figure D-iO, a histogram is presented of TNT equivalencies based on the 
URS terminal yield data tabulated in Appendix A. The mean value of 
the TNT equivalencies for these data is 20.4%. It should be noted 
that these data show that there is a large concentration of very low yields. 

This contradicts one of the criteria specified for the density function for the 
previous cases; that is, that values of F(z) in the vicinity of zero yield 
are relatively low in comparison with values at higher yields. A probability 
density function was then developed based on these data for use in the analysis. 

It should be noted that no effort was made in this portion of the study to verify 
that these test data actually represent a good statistical sam* 'e for application 
to this type of analysis. For instance, the experimental data were obtained 
primarily from tests with very small quantities of propellants and the tests 
were run under conditions which will not exist in the actual vehicle. For this 
reason, the reader is cautioned against placing any special significance to the 
results using this model. 
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Figure D-8. Composite Probability Density Function - Model 3 





Figure D-9. Probability of Exceeding a Given Yield - Model 3 








The Beta function was also used to define the density function for this model. 

In this case the equation is 

F(z) = 6. 5 (1 - z) 5 * 5 

1 

The density function defined by this equation has a mean value of 20% as was j 
noted for the experimental data. The density function is plotted in Figure D-ll. 
The correlation between the assumed density function and the experimental 
data is shown in Figure D-12. In this figure, the horizontal scale is in % TNT 
equivalency with a maximum value of 150% as was used in Figure D-10 for 
the experimental data. It should be noted that the F(z) values for the assumed 
density function are modified in this figure so that the ?rea under the curve 
remains unity. 

In Figures D-13 and D-14, the results of the analysis using the distribution 
model shown in Figure D-ll are presented. The results of the analysis, show 
that the mean value of yield would be 0. 19 X 10^ lb, the 95% value would be 
0. 94 X 10^ lb and 0. 90 X 10^ lb for Model 3 for which the probability density 
function was assumed without benefit of the test data. However, the mean 
values of the density functions in both models were the same; i. e. , 20%. 

The experimental data summarized in Figure D-10 have a mean value of TNT 
equivalency of approximately 20% and a maximum value of approximately 100%. 
However, other test data and theoretical studies indicate that even higher 
equivalencies may occur. Values as high as 300% are considered possible by 
some sources. 

As previously discussed, a maximum value of the TNT equivalency of 100% 
was first assumed (Model 1). Subsequent calculations (Models 2, 3, 4) were 
based on a maximum value of 150%. The sensitivity of the Beta distribution 
function, which was developed based on the experimental data (Model 4), is 
now summarized for changes in the assumed maximum value of yield. One of 
the features of the Beta function is that it Las, in each case, an upper limit 
(i. e. , it does not admit to values exceeding this limit). The other general 
characteristics of the Beta distribution function shown in Figure D-ll were 
retained, including the mean value of 20%. 
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Figure D-12. Comparison of Assumed Probability Density Function, 
Model 4, to Experimental Data 
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Figure D-13. Composite Probability Density Function - Model 4 




Figure D-14. Probability of Exceeding a Given Yield - Model 4 



In Figure D-15, the density function for Model 4 from Figure D-ll (based on 
150% maximum TNT equivalency) is shown together with the function for the 
cases where the maximum value of the distribution is increased to 300% and 
infinity. As indicated, there is no discernible difference between the 300% 
case and the infinity case and only minor differences between these cases and 
the one for a 150% limit. The insignificance of the differences in the data 
presented in Figure D-15 is further illustrated by comparing the yield/ 
propellant weight values for certain cumulative values of probability shown 
below: 


Yield /Propellant Weight Value for 
Maximum Yield Indicated Cumulative Probability 



80% 

90% 

95% 

1. 5 X Propellant Weight 

0. 329 

0.448 

0.554 

2. 0 X Propellant Weight 

0.328 

0.451 

0. 566 

2. 5 X Propellant Weight 

0. 327 

0.454 

0. 574 

3. 0 X Propellant Weight 

0. 326 

0.455 

0.578 

00 

0. 322 

0.461 

0. 599 


The above data indicate that the yield /propellant weight value corresponding to 
a 95% cumulative probability changes approximately 4% when the maximum 
value (y n ) is doubled and increases approximately 8% for the infinite upper 
yield case. Thus it is concluded that when the mean value of the density 
function is established, the effect of the assumption of larger maximum values 
of yield will be negligible when a Beta function with the general characteristics 
shown in Figure D-ll is used. 

D. 3. 5 Model 5 

In previous sections, a Beta function was used to define the probability that 
a specific yield would result in event of an explosion. A basic characteristic 
of this function as used in the preceding analyses is that the probability of 
occurrence goes to zero at the maximum yield value. 
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Figure D-15. Comparison of F(z) Values 




The question can be asked as to the effect of a distribution which v/ould assume 
some finite probability that the maximum value will be attained. To evaluate 
this effect, a uniform distribution of the probability as a function of yield 
(Figure D-16) was assumed together with a maximum value of 150% TNT 
equivalency. Under this assumption, there is the same probability (for a 
given Az value) that the yield in event of an explosion will be zero, maximum, 
or some value in between. This distribution is considered very unrealistic 
when one considers the more or less optimum conditions relative to mixing 
etc. that must exist to obtain high yields. It is therefore considered as an 
extreme limit on the range of probability distributions that might be considered. 
In Figure D-i7 the composite probability density function is shown v<;sed on this 
assumed distribution on yield • t an explosion and the probability o. various 
amounts of propellar.t being involved in the explosion. The expected ralue for 
this distribution occurs at yield of 0. 67 X 10^ lb with a 95% cumulative 
probability occurring at approximately 3. 5 X 10^ lb. 

In Figure D-18, the cumulative probability of thcyield exceeding specific values 
is presented as a function of yeild. This figure show* that the probability of 
exceeding a yield of 1 million lb unoer this model is 15%. As would be expected 
because of the characteristics of the distribution assumed, this value is higher 
than that for the other assumed models. 

It should be noted that with the uniform probability model there is approximately 
one chance in five of exceeding a yield of 0. 8 X 10^ lb (20% of the 4 million lb 
of propellants). This value is similar to that for Models 1 and 2 in which a Beta 
function was used which specifies zero probability that a zero yield or a 
maximum yield will result. The above comparison shows that the most signifi- 
cant difference ia in the probability of attaining high yields. This results from 
the fact that the uniform distribution achnowledges a significant probability that 
maximum yields will occur. It should also be kept in mind in reviewing the 
above results that the uniform distribution model is presented for comparative 
purposes as a limiting case. 
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Figure D-17. Composite Probability Density Function - Model 5 
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Figure D-18. Probability of Exceeding a Given Yield - Model 5 


\ 4 



D. 4 SUMMARY 

In Table D-2, the results of the analysis are summarized. It shows that for any 
of the density functions assumed there is a relatively low probability (conditional) 
of attaining yields approaching those corresponding to 60% TNT equivalency 
based on all the propellants aboard the space shuttle. The total probabilities of 
attaining the indicated yields would be many orders of magnitude less. 

It is interesting to note that the lowest probabilities are obtained using Models 3 
and 4. Model 3 was based on a Beta function which had zero probability of 
attaining zero yield or maximum yield and had a mean value of 20%. Model 4 
was based on the experimental data which showed a high probability of attaining 
low yields. 


The highest probabilities are generally associated with the uniform distribution. 
It is emphasized that this is not considered to be a realistic model since it 
assumes that there is the came probability of attaining any yield including zero 
and the maximum value. 


Table D-2. Comparison of Probability Models Investigated 


Probability of Exceedin 



Model 

0. 8 X 10° lb 

2. 4 x 10 b lb 

4. 0 X 10 6 lb 

1 

y max = 100 % : 7 = 60 % 

0. 19 

0. 04 

0 

2 

y m ax = 150%; V = 60% 

0. 19 

0.06 

<0. 01 

3 

y max = 150%; y = 20% 

0. 06 

<0. 01 

<0. 01 

4 

(Exp. Data; y max - 150%) 

0. 06 

<0. 01 

<0. 01 

5 

(Uniform Distribution 

0. 19 

0.09 

0. 03 


y 


max 


150%) 
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